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1 Introduction and overview 



One of the most exciting discoveries in the last few years is the integrability of the 
maximally supersymmetric Yang-Mills theory [H EJ [3] and its relation to the superstring 
theory in AdS 5 x S 5 background |U El IB]- An all-orders version of the Bethe Ansatz 
equations for AdS/CFT, valid asymptotically, was first formulated by Beisert, Eden and 
Staudacher [7J, based on previous work [HI El CD]- The proposal by Beisert, Eden and 
Staudacher solves the crossing symmetry conditions, formulated by Janik [Uj. This 
proposal, interpolating all the way from strong to weak coupling, was intensively tested, 
the most sophisticated tests being performed on the large spin limit of the so-called twist- 
two, or cusp, anomalous dimension. In the limit of large Lorentz spin S, this quantity 
scales logarithmically [HI H31 HI HS1 CB] 

A-S = f(g)\nS + ... , (1) 

where g is the coupling constant, related to the 't Hooft coupling constant A by 

The universal scaling function f(g) was computed perturbatively in the gauge theory 
up to the fourth order. The third order result was extracted [UJ from a QCD computation 
by Moch, Vermaseren and Vogt [18J. The universal scaling function appears in the 
iterative structure of the multigluon amplitude [TH] and it was computed to the third 
odrer in [20] and numerically to the fourth order, after an impressive effort [21, 22J. 
On the string side, the universal scaling function was also computed for the first three 
non-trivial orders [231 |23| |25] [26] 

ft \ a 31 °g 2 K 1 (>>\ 
f(g) = ±g -—- + ..., (3) 

7r Att z g 

where K= (3(2) is Catalan's constant. It is remarkable that both the weak coupling and 
the strong coupling results for the universal scaling function can be reproduced from the 
conjectured Bethe ansatz equations. In this context, it is determined by the integral 
equation, written down by Eden and Staudacher [T5] 



1 f°° a(v) f°° ( 1 



"(«) = - J _ ^ {u _ v)2 + 1 ~ J _ dv K(u, v) [^)-—). (4) 

With the integration kernel K(u,v) determined in [7J, this equation is known as the 
Beisert, Eden and Staudacher (BES) equation. The universal scaling function is given 
by the integral of the density [27] 



f(g) = 16g 2 J a(u)du. (5) 

Although it is not, at least for the moment, possible to solve the equation (j3J) in 
a closed form for arbitrary g, it is relatively easy to extract from it the perturbative 
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expansion at weak coupling [T5l [7] . This coefficients of the perturbative expansion agree 
with the field-theoretic results [17] 1201 [2T| 122] . 

The strong coupling limit of the equation (jl]) proved to be much more difficult to 
master analytically. The first results were obtained numerically [2H] arid they correspond 
to the three coefficients in ([3]). The first coefficient in ([3]) was obtained analytically by 
various methods (23, [29j [30j EI] , while the second was obtained, although not from the 
BES equation, by Casteill and Kristjansen [32] and later by Belitsky [33]. The reason 
the strong coupling limit of the BES equation is so difficult to take is that the expansion 
of the scattering phase in powers of 1/g is not uniform in the rapidity variable u. There 
are three different regimes for u which are to be considered. The first is the plane-wave 
limit [31], where \u/2g\ > 1, or in terms of momenta p ~ 1/g. The second is the so- 
called giant magnon regime [35], with \u/2g\ < 1 or p ~ 1. The third regime was called 
[3"o] the near-flat space regime, at it correspond to to u ± 2g ~ 1, or to momenta of 
the order p ~ 1/ \fg~- As Maldacena and Swanson pointed out in |36j, in this region 
the perturbative expansion of the dressing phase is completely reorganized, compared to 
that in plane-wave and giant magnon regimes. This is the reason why the attempts to 
to solve the BES equation order by order in 1/g failed beyond the leading ordei0 . 

A big step forward was made recently by Basso, Korchemsky and Kotahski [37], who 
succeeded to give a procedure for obtaining all orders in the strong coupling expansion 
recursively. One of the important ingredients of their work was to linearize the Bethe 
ansatz equations by transforming the BES equation into a set of two equations, by 
exploiting the expression of the dressing kernel as a convolution of two "undressed" 
kernels [7J. The result is a set of two coupled integral equation, for the physical and an 
auxiliary density. The idea of linearization was first proposed by Kotikov and Lipatov [27] 
and subsequently by Eden [38J. Basso, Korchemsky and Kotahski [37] used the Fourier 
representation of the BES equation, which is more adapted for numerical analysis, as 
well number of numerically inspired hypotheses. 

The present work arose as an attempt to derive the results of [37] by purely analytical 
consideration. Our previous work [30J on the strong coupling limit of the BES equation 
suggests that a natural way of solving this equation is to reformulate it as a functional 
equation for the resolvent. Here, we will pursue this direction. To use the resolvent 
to solve the Bethe equations is by no means a new idea; in particular it was used in 
[32| 133] . The novelty of our approach is to replace the density in the BES equation (TJ]) 
by the resolvent and to exploit the analyticity property of the latter. After writing the 
linearized BES equations in terms of the resolvent, it is possible to transform the integral 
BES equations into a set of functional equations for the physical and auxiliary resolvents. 
This procedure works for an arbitrary value of the coupling constant]^ 

The question if the functional BES equations can be solved for any value of the cou- 
pling constant is still open. What is clearly possible to do is to solve these equations 
at strong coupling in perturbation theory, that is, by neglecting the non-perturbative 
correction. This is possible because, in the absence of non-perturbative terms, the resol- 

the approaches |32[ 133]. which in principle treat a different limit than the BES equation, the 
contribution from the near-flat space regime is suppressed from the beginning. This rapidity regime 
falls inside the gap of the density. 

2 The essential steps of this procedure are already present in the paper [57] by Kotikov and Lipatov. 
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vents posses extra symmetry properties. The algorithm of solving the equations loosely 
follows the one of Basso, Korchemsky and Kotahski [37J. 

First, we find the general solution for the resolvents at strong coupling in the giant 
magnons/plane-wave regimes. In these regimes the rescaled rapidity u/2g is kept finite. 
The functional equations are linear and homogeneous, so that the general solution is a 
linear combination of a countable set of particular solutions. Each of these functions 
have a non-integrable singularity at the points u/2g = ±1. The indeterminacy in the 
coefficients of the linear combination and the singularities of the individual solutions can 
be cured by analyzing the solution in the near-flat space region, where the perturbation 
series in 1/g is reorganized compared to the giant magnon and plane- wave regimes. 
In the near-flat space regime the resolvents must have a series of integrable square root 
singularities in the lower half-plane and must decrease as 1/u at infinity in the upper half- 
plane. Remarkably, the condition that the solutions in the plane- wave/giant magnon and 
near-flat space regimes match analytically allows to determine uniquely the resolvents in 
the strong coupling limit up to non-perturbative corrections. This matching condition 
is equivalent to the quantization condition by Basso, Korchemsky and Kotahski [37]. 

We formulate a recursive procedure for computing analytically the coefficients of the 
1/g expansion of the density of Bethe roots. In the plane-wave regime, \u/2g\ > 1, the 
density has a fourth order branch point at u/2g = ±1 at any order in 1/g. 

The paper is organized as follows: in section 2 we derive the linearized equations for 
the resolvents, in section 3 we transform the integral equations into functional equations 
and we discuss the analyticity properties of the resolvents and in section 4 we find 
the perturbative solution by imposing the required analyticity properties to the general 
solution. 



1.1 Notations 

The notations used in this paper are similar with those of our previous paper |30j. We 
will denote by e the inverse gauge coupling and will use a rescaled rapidity, more adapted 
to the strong coupling limit: 

1 ^old /„n 

as well as the variable x(u), related to u by 



Note the branch cut of x(u) for u £ [—1, 1] . In the intermediate steps we will also use 
the notations 

x ± {u) = x(u ± it) . (8) 

Sometimes it will be useful to switch to the parametrization which resolves the square 
root of x{u) and which is the hyperbolic limit of the elliptic parametrization in [30J: 

u = coths, x(s) = coth-, = e 2s . (9) 
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The BES equation is formulated for the density function u{u) in the limit of large spin 
S 3> 1, related to the distribution p{u) of Bethe roots by 

p{u) = (^-a(u)^j logS + O(S ), u<^S/2g. (10) 

2 Linearized BES equations for the resolvent 

In this section, we reformulate the Beisert, Eden and Staudacher equation [USE], which 
determines the density of Bethe roots corresponding to the twist-two operator, as a set of 
two equations for the physical resolvent and an auxiliary resolvent. This is essentially the 
program which was carried out by Kotikov and Lipatov [2Z] and by Eden [3H] , although 
they have not explicitly identified the inverse Fourier transform of the density on the 
positive half- axis as the resolvent. 



2.1 The resolvent 

Let us consider the density a(u), supported on the real well as its Fourier 

transforrr|§ a(t) 

/oo 
due itu a{u). (11) 
-oo 

The resolvent is defined, as usually, by 

oo 

iWO)= / dv^-. (12) 
J u - V 

— oo 

and it can be seen as the inverse Fourier transform of the density a(t) on the positive 
half-axis. Let us assume that u is in the upper half plane. Then we can write 

/oo poo 
dv / dte it{u ~ v) a(v) (13) 
■oo JO 



oo 

itu 



dt e ltu a(-t) . (14) 
The density is given by the discontinuity of the resolvent across the real axis: 

(?(u) = [R phys {u - iO) - i? p h ys (w + i0)] (15) 

Since the density is supported by the whole real axis, the resolvent is given by two dif- 
ferent analytic functions in the upper and in the lower half-planes. Due to the symmetry 
a{— u) = o(u) we have the relation 

-Rphys(-W) = --RphysH (16) 



J Our definition of a(t) is slightly different of that of [T5] and the subsequent references. 
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Therefore, the symmetry determines the resolvent in the lower half plane once the resol- 
vent in the upper half plane is known. The behavior of the resolvent at infinity is related 
to the universal scaling function 

oo 

(17) 



2.2 Linearizing the BES equation 

The linearization of the BES equation is best performed on the Fourier transformed form, 
although the manipulations can be done abstractly without reference to a particular 
representation. In our notations, the Fourier transformed BES equation reads 

r°° dt' 

(1 - e-™)a{t) = - — {Kit, t') + K d (t, t')) (a(t f ) - a (t')) . (18) 
Jo 27r 

where <Jo(t) = 4eS(t). Here we use the conventions 

K{t,t') = 2nte-^J2 2n ^¥^- = K+{t,lf) + K^tJ) , (19) 

where K + and K_ contain the expansion on odd and even order Bessel functions respec- 
tively, and the dressing kernel is given by the "magic formula" [7] 

f°° dt" 1 
K d (t,t>) = 2 —K-(t 1 f) 1 _ e _ Mf , K + (f 1 lf) • (20) 

Two representations of the dressing kernel in the rapidity space were given in [30J and 
[39] . It is convenient to introduce an operator SS, diagonal in Fourier representation: 

= 1^=1* • (21) 

In the rapidity space 

3S~ 1 = 1-D, where Df(u) = f(u + 2ie). (22) 

The equation ( fl8|) can be written symbolically as 

- 2a = [(1 + 2SSKJ)(1 + 2SSK+) - l](a - a ) . (23) 

Now it is possible to transform the BES equation into a pair of equations with the 
"main" kernels K± appearing linearly. This can be done at the expense of introducing 
an auxiliary density r defined by 

T + a = -(l + 2SSK + )(a-a ). (24) 

The linear system of equations obeyed by the physical and auxiliary density a and r is 
simply 

r + a = -2SSK + (a - a ) (25) 
t-o = -2SSK_(r + a ) , 
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00 dt' 

00 dt' 
2~ix 



K + (t,t') (a-a )(t') 

K_(t,t')T(f). 



(26) 



or, in Fourier representation, 

(l-e-™)(r(t) + a(t)) = -2 

(1 - e-^)(r(t) - a(t)) = -2 
In the last line we have used that K_(t, 0) = 0. 

2.3 Holomorphic BES kernels 

The next step is to transform back the equations (I2"6"j) in the rapidity space. As we men- 
tioned above, the inverse half-space Fourier transform of the density a(t) gives the resol- 



vent -Rphys(w)- The latter defines a pair of functions R up (u) and R n (u) = —R up (—u), 
analytic respectively in the upper and lower rapidity half-planes. Because of the sym- 
metry property, we are going in the following to concentrate exclusively on R np (u). It is 
this function, together with its analytical continuation beyond the real axis, which will 
be denoted in the following by _R p hy S (w). 

Assume that > and perform the half-space inverse Fourier transformation to 
the rapidity plane, 



00 00 

dt e ^ [ 5k( M ')/(*') 



2tt 



2vr 



dv 



dt 
— < 
2vr 



:it/u 



dt 
— < 

2vr 



-it'v 



K(t,t') 



f>7(0 • (27) 











Therefore, since we intend to work only with functions defined in the upper half plane, 
we can retain only half of the original kernel in rapidity space, namely 



K(t,t'). 



(28) 



Here we use the superscript e for the kernel in order to indicate that it depends on the 
coupling constant g = l/4e. Explicitly the "holomorphic" part of the odd and the even 
kernels reads 



K e (u,v) 



K%{u,v) 



1 d 


In 


'1- ' 


2m du 


V x+ y 


1 d 


In 


\- 1 


2iri du 





+ In 1 + 



-In 1 + 



x + y~ 



(29) 



The dependence on e in ( !29l) comes only through the shifts x ± = x(u ± it) and it will be 
removed by change of variable and shift of the integration contour. The e — > limit of 
the kernels ( |29l) will be denoted without superscript 



K~(u, v) 
KAu,v) 



2i\i 1 — x" 1 



1 



1 



2m' 1 — x 2 



y + x 
1 



y + 



(30) 
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When using these equations, one has to remember that > and ^sv < 0, or 

x = x(u + i0), y = y(v — iO) — l/y(v + iO) . (31) 

2.4 BES equations for the resolvents 

For later convenience, we introduce the shifted resolvents 

oo 


oo 

H(u) = -i I dt e iut e et r(t) . (32) 



o 



This definition is valid for Q[u] > and the first singularity for R(u) and H(u) is situated 
on the real axis. R and H can be analytically continued to Q[u] < 0. We will introduce 
another pair of functions by 

R±(u) = l[R(u)±H(u)\ (33) 
as well as the related functions r±(u) 

r± (u) = R± (u) -R±{u + 2ie) . (34) 



Now we can take the inverse Fourier transform of the equations (1261) and make the 
shifts u — > u — ie and v — > v + ie, in order to get rid of the e-dependence in the kernels 
K± (u,v). We obtain the equations^ 



Aie 

R+(u) - R+{u + 2ie) = 



x 2 



- - J dv K+(u,v) [R+{v + 2ie) + R.(v + 2ie)\ 

R„{u) -R_{u + 2ie) = J dv K„(u, v) [R + {v + 2ie) - R_{v + 2ie)} . (35) 

After the change of variables the integration contour for v goes along the shifted real 
axis M. — ie, but it can be placed anywhere between the branch cut [—1, 1] of the kernels 
K± and the branch cut [—1 — 2ie, 1 — 2ie] of the resolvent R(u). We will assume that the 
contour is just below the real axis. The variable u originally lies in > 0, but we can 
analytically continue it to the whole rapidity plane, using that the integration kernels 
are holomorphic. 



3 Functional equation 

The kernels (I3TJ1) look almost like Cauchy kernels, if not for the branch cut of the variables 
x(u) and y(v). This suggests that we may simplify the BES equation further. This can 
be done provided that we know the analytical properties of the functions the kernels 

This way of rewriting the BES equation have been first proposed by Kotikov and Lipatov in [27) . 
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act on. In this section, we derive the analytic properties of the resolvents R±(u) and of 
the functions r±(u) and translate the action of the kernel in terms of an integral on the 
interval [—1,1]. This transformation will allow to transform the integral equations into 
a functional equation. 

3.1 Analytic properties of the resolvents 

We start with the linearized BES equations 

r+(u) = -j— - - / dvK+(u, v) [R+(v + 2ie) + R.(v + 2ie)] 

R—iO 



riu) 



where 



J dvK_(u, v) [R + (v + 2ie) - R-{v + 2ie)\ (36) 



/ s 12/1 1 . 

K±(u,v) = t r± r 37 

2ni 1 -x 2 \y- i y + - 1 



with x = x(u) and y = x(v). The variables u and v belong to the physical sheet, which 
means that \x\ > 1 and \y\ > 1. Since the kernel becomes singular only for u in the 
interval [—1,1], and there is no other singularity when u and v are on the physical sheet, 
the functions r±(u) are analytic in C\[— 1, 1]. We deduce that the resolvents 

oo 

R ± (u) = r ±{u + 2ine) , (38) 

n=0 

have a semi- infinite set of equidistant cuts as shown in Fig. 1. From the explicit form of 
the kernels it follows that 

r + (u) oc , r(u) oc 4r (u — > oo) . (39) 
u u 

The large u behavior of the resolvents is R+{u) ocl/u and R_(u) oc 1/u 2 . 



-1 




u 


1 


-1 


u 

1 










lie J 































Fig. 1: Left: Physical sheet for r±(u). Right: Physical sheet for i?±(u). 
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3.2 Analytic properties of the kernels. 



Since the resolvents have definite analyticity properties, we can considerably simplify the 
action of the integration kernels. Consider the integral 

+00— tO 

J dv K±(u,v)F(v), (40) 

— 00— iO 

where F(v) stands either for R + (v + 2ie) or for R-(v+2ie). The function F(u) is analytic 
in the half-plane Qu > — e. We will actually need a weaker assumption, namely that the 
function F(u) is analytic in the upper half plane > with the real axes included. 

We would like to place the integration contour above the real axis, since in the upper 
half-plane both K± and F are analytic. Using the properties 

x(v-iO) = l/x(v + i0) , u e [-1,1] 

x(v-i0) = x(v + iO) , u G K\[— 1, 1] (41) 

we can place the integration contour above the real axis at the price of changing the form 
of the kernel in the interval [—1, 1]. Next, since in the upper half-plane the function F 
is analytic, we can deform the rest of the contour to to a contour that goes along the 
interval [—1,1] in the opposite direction. Adding up the two contributions we evaluate 
the integral as 

1+iO 

Tr ril . 2 f dv ri . . / — yx yx 1 1 \ 

K±F(u) = / F(v) ± — r =F r 

KJ 1-x 2 J 2m KJ \y-x y + x y-l^y + lj 

-l-KQ 
1+tO 1 

/dV V ~ y f 1 1 \ 

2^i x- ±\v-u T v + u) ' ^ ' 

-1+iO x 

In the original equation (|35l) . the integration contour is pinched between two cuts dis- 
tanced by e and the limit e — > is not well defined. On the contary, the form f)42p of the 
integration kernel has a smooth e — > limit: 



dv ll-v 2 F(v + iO) ± + iO) 
2n\u 2 — l v — u 



K ± F = /^\/^-r- r— f ( 4 3) 



-1 



3.2.1 Projective properties of the kernels K± 



Via the transformation v —v the integral in (1431) can be written as an integral over a 
closed contour around the segment [—1, 1]: 



— F(v) 

2m' V u 2 — 1 v — u 



K ± F = <P — F(v) \l 7T— 7 — (44) 
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with 




if > , 

F (u) = { , J \ . fCV _ (45) 

it < . 

Denote by £ ± the linear space of even/odd functions, analytic outside the interval 
[—1,1] and decreasing at infinity faster than 1/u at infinity. For any f± e , the kernel 
iT-t- acts as the identity operator: 

K±f± = f±, f±eC ± . (46) 

In particular, 

K + r + = r+, i^_r_ = r_ . (47) 

Now consider the result of the action of the kernel K± on an arbitrary function F(u). 
Since the kernel K±(u,v) is even/odd as a function of u, is analytic outside the interval 
[— 1, 1], and decreases as 1/u 2 at infinity, the resulting function F± = K±F belongs to 
C ± . As a consequence, K\F = K±F± — F± — K±F. We conclude that the kernel K± 
is idempotent: 

Kl = K ± . (48) 

An example of a function that does not belong to is the constant function. The 
action of K±, evaluated by expanding the contour to infinity, is 

K+-l = —^= = -^, K.-1 = 0. (49) 
V« 2 - 1 1- x z 



3.3 The BKK transformation 

The linearized BES equations (135]) considerably simplify when written in terms of the 
functions r + and T_ defined bjo 

r+(n) + r_(n) = R + (u) + R-{u + 2ie) + 2ie 

r+(u)-r_(u) = R-(u) -R+{u + 2ie) + 2ie. (50) 

Indeed, with the help of the identities (j4"9l and ( 1471) we write the linearized BES equations 
(1551) as 

ir + (r + + r^) = o, (5i) 
Ar_(r+ -r_) = o. (52) 

Therefore the solution of the BES equation is a linear combination of the zero modes of 
the operators K±. 

s This transformation corresponds to the one given by eq. (6) in [37] . 
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We have seen that the action of the kernels K± on the functions analytic in the upper 
half plane is given by (1431) . As a consequence, the necessary and sufficient condition that 
the function F is annihilated by K± is 

F(u + iO) ± F(-u + iO) = , ue [-1, 1] . (53) 

Note that the condition (j53p does not imply the function F is odd or even. It can be 
written in terms of the variable 

F(x) = TF{-l/x) . (54) 

For example the first equation (1491) means that fo(x) = is a zero mode of K + . This 
function satisfies /o(— 1/^) = ~fo( x )- 

3.4 From integral to functional equations 

Now we can reformulate the homogeneous integral equations flSTj) as a pair of functional 
equations for T + and T_. According to ( 1531 or ( 1541) . these equations imply the following 
boundary conditions on the upper edge of the cut [—1,1], 

T + (u + iO) + T_(-M + iO) = 0, we [-1,1], (55) 

or, in terms of the variable x, 

r+(-i/x) = -r_(x). (56) 

The last equation should hold on the arc \x\ = 1, Qx > 0. 

Hence the solution of the BES equation must be among the solutions of the functional 
relation. Note that this equation is exact in the sense that in the derivation we did not 
assumed that e is small. 

Of course this relation has a huge set of solutions. The physical solution is dis- 
tinguished by imposing its analyticity properties in the vicinity of the singular points 
u — *■ oo and u — ±1. The extra conditions that single out the physical solution are 
formulated in terms of the original resolvents R± = \{R ± H), or equally the functions 
r±{u) defined by (|3~4"|) . Namely, the functions r± must be analytic everywhere outside the 
cut [—1,1], where they have square-root singularities, and behave at infinity according 
to ( |39l) . This conditions determine the analytic properties of T±, which are related to r± 
by 

oo 

T_(u) = \r + (u) — |r_(w) + r + (u + line) , (57) 

n=l 
oo 

T+(u) = lit + \r + {u) + |r_(u) + J^r_(w + line) . (58) 

n=l 

Since the functions r± have a square root cut along the interval [—1, 1] of the physical 
sheet, all singularities of the functions T± are of square root type. 
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4 Perturbative solution at strong coupling 



In this section we obtain the perturbative solution for the resolvent. First we consider 
the limit e — > with u fixed. This limit corresponds to the plane waves (PW) or giant 
magnons (GM) regimes, depending on the interval where the rapidity takes its values 
(Fig. 2). The distribution of Bethe roots in the PW and the GM regimes is given by two 
different analytical expressions, but for the resolvent they they are related by analytical 
continuation. 




Fig. 2: The physical density p(u) = 2e/n — cr(u) in the strong coupling limit and the three 
regimes: plane waves (PW) for u < —1 and u > 1, giant magnons (GM) for —1 < u < 1, and 
near flat space (NFS) in the vicinity of the points u = ±1. 



It happens that in the strong coupling limit, and in all orders in e, the intricate 
cut structure of the resolvent and the related functions can be replaced by a single 
cut u G [—1,1], but with fourth order instead of second order branch points at u — ±1. 
Furthermore, an important simplification stems from the fact that in the PW/GM regime 
the combinations T±(u) have definite parity, 

r±(-u) = ±r±(u) . (59) 

This will allow us to write the general solution of the linearized BES equation. 

Since the equations are homogeneous, the general solution is a linear combination of 
all particular solutions with arbitrary coefficients c n , which are functions of the coupling 
constant. The behavior of the resolvent at u —>■ oo gives one linear constraint on the 
coefficient functions c n (e), which is not sufficient to determine them. 

The rest of the information is supplied by the conditions on the analytic propertirs of 
the solution in the vicinity of the singular points u = ±1. For this purpose we blow up the 
vicinity of the the two singular points so that the cut structure of the resolvent reappears. 
Instead of keeping u fixed, we take the limit e — > either with z — {u — l)/e fixed or 
with z = (u — l)/e fixed. This strong coupling limit corresponds to the near flat space 
(NFS) regime [36]. Then we compare the power series expansion at z = 0, which follows 
from the analytic structure of the exact solution, with the expansion at z = oo, which 
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is determined by the perturbative solution in the PW/GM regime. The requirement 
that the two expansions match with each other is sufficient to determine both of them, 
order by order in e. Technically it is more advantageous to compare the inverse Laplace 
transforms for which the shift operator D becomes diagonal. A recurrence procedure, 
analogous to that of [37] , allows to obtain analytically the density of Bethe roots in any 
order in e, both in the PW/GM and NFS regimes. We check that the result of [37] for 
the universal scaling function is correctly reproduced. 

4.1 General form of the solution in the PW/GM regime 

Let us first prove the symmetry property ( 1591) . For that we use the expression of T± in 
terms of the functions of definite parity r±{u) = ±r±(— u), given by (1571) and (1581) . We 
observe that the combinations 

I±(u) T r±(-u) = T^r T (« + 2me) (60) 

are periodic functions with period 2ie. From here and from the fact that r±(u) ~ 1/u 2 at 
infinity it follows that the r.h.s. (|6"U1) vanishes in the limit e — ► up to non-perturbative 
terms. To see that we perform Poisson resummation. Assuming that > 1, we have 



V r±(u + 2ine) = - V e~™ u/e <b ^L e mv/e r ± (v) , 
— ' e — ' / 2m 

=— oo n=l 



(61) 



where the integration contour closes around the physical cut [—1,1] of r±. The series in 
g-vrnu/e j g ra pj(jiy convergent when u > 1 and diverges at u = 1. When Jiw < — 1 we get 
a similar expansion, but the opposite sign in the exponents. In both cases the result is 
exponentially small except at the points u = ±1. 

Therefore, if we neglect these non-perturbative corrections, the solution should have 
the additional symmetry (1591) . Then the functional equation (1551) can be replaced by a 
simpler one, 

r + (u + io) = +r_{u-io) 

r+(u-io) = -r_(u + io). (62) 

We remind that these equations are valid on the cut, where u G [—1, 1]. From here it 
follows that r± are obtained as different branches of the same meromorphic function, 
defined on a four-sheet Riemann surface. 

In terms of the global parameter of the Riemann surface, s, the functional equations 
get the form of periodical conditions 

T±(s±in) = ±r T (s). (63) 

It is convenient to work with the combinations 

g± = r + ± *r_ , (64) 
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for which ( |59i) and ( |63i) take the form 

G±(-s) = G T (s), G ± (s + ztt) = ±iG±(s). (65) 
We can represent the general solution in the form of the series: 

G±(s) = 22e^c„(e)e ± ( 2n+1 / 2 ) s 

^1 • (66) 

We will first obtain the general form of the solution in the PW/GM regime the branch 
points condense into continuous lines starting at the points u = ±1 and the resolvents 
are described, as we will see later, by meromorphic functions with a single pair of branch 
points at u = ±1. 

The perturbative solution fl66|) is valid in the limit where the distance between the 
subsequent branch points vanishes and the infinite sequence of simple branch points 
starting at u — ±1 produces a fourth order branch singularity at u — ±1. 

The solution has three singular points, u = l,u = — 1 and w — > oo. As usual in such 
kind of problems, the coefficients functions c n (e) in the series ( 166]) will be evaluated by 
matching with the asymptotic behavior of the solution at the singular points. 



4.1.1 The density in the GM regime 

The general form of the solution given by (IfiT)]) is sufficient to determine perturbatively 
the density in the giant magnon regime. Indeed, inspecting each of the terms, one can 
verify that the value of the resolvent on the interval — 1 < u < 1, and therefore the 
density, is constant and is given by the leading order. 

This fact is actually a direct consequence of the equations ( 1631) and the symmetry ( 1591 . 
Indeed using the anti-symmetry of the resolvent -R p h ys , we can express the fluctuation 
density in terms of the values of the resolvent above the real axis: 

a(u) = -^-[R phys {u + iO) + R phys {-u + iO)} (67) 

= -7^[R(u + ie) + R(-u + ie)}. (68) 
Further, by the definition ( |50l) . the resolvent -R p h ys is expressed in terms of T± as 

R phys = -2it + - 2 £) _ 1 (z^r_ + zHr + ) (%u>o) 

where D is the shift operator defined by (!22l) . Applying the functional equation ( J62l) . we 
see that all the terms on the r.h.s. of ( IBTj) except the constant term cancel and therefore 
to all orders in e 

a(w) = 2e/7r, we [-1,1]. (69) 

This means the distribution of Bethe roots has a gap on the interval [—1,1]. The physical 
density ( jTUi) . which gives the distribution of Bethe roots, vanishes to all orders in e in 
the GM regime. 
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4.1.2 Expansion at u = ±1 and a scaling condition for the coefficients 



Let us examine the behavior of the solution (1661) near the singular points u = ±1. 
We mentioned that the strong coupling limit is not uniform in u. The strong coupling 
solution have different properties in the limit considered above, 



and the limit 



e^O with u fixed (PW/GM), (70) 



u 2 - 1 

with fixed (NFS) . (71) 



The singular behavior at u — ±1 in the PW/GM limit is an artifact of the rescaled 
rapidity ([6]). If we take the NFS limit ( 17TT) . the solution for the density must be integrable 
at it = il. It is obvious that the strong coupling expansions in the two limits do not 
match since the solution (|66|) gives non-integer powers of e when considered near u = il. 

Our analysis of the analytical properties of the solution allows us to determine its 
general form near u = il. The conditions that it goes smoothly into the solution (1BT)|) 
obtained for the rest of the complex plane will be used in the next section to fix the 
coefficients c n . 

The complex variables relevant for the vicinity of the points u — 1 and u = — 1 are 

u — 1 u + 1 

-~ = — < 72 > 

The variable z coincides, up to a shift by 2g, with the original (before rescaling by 2g) 
rapidity in the BES equations. 



z 





II 





















Fig. 3: Physical sheet for R±(u) in the NFS limit. 



In the NFS limit the cuts become semi-infinite, with the branchpoints placed at at 
z = 0, — i, —2i, . . . , as shown in Fig. 3. The functions r±(z) have by construction an 
integrable square root singularity at z — 0. Therefore they can be expanded at small z 

as 

r±(z) = b n& * n ~ 1/2 + E zn ( W < 1 )• (73) 

n>0 n>0 
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The compatibility of the expansions ( 1731) and ( |66l) imposes strong restrictions on the 
coefficient functions c n (e). In particular, each term of the expansion (166]) must have a 
non-singular limit e — > when expressed in terms of the variable z or z. This means 
that the coefficients c n (e) must scale as e' ra ', so that their Taylor series have the form 

oo 

c„(e) = e |n| a„(e), a n (e) = oc n , p e p (n 6 Z) . (74) 

p=0 

We arrive at the following expression of the general solution in terms of the shifted 
rapidity variable z: 

— -) ■ (75) 

The strong coupling expansion of G± with z kept fixed is different than the expansion 
of the solution with u fixed, (1661) . In particular, it contains fractional powers of e. The 
resolution of this paradox is in the non-uniformity of the strong coupling expansion with 
respect to the rapidity variable u. Near the singular points u = ±1 the strong coupling 
expansion should be performed according to the prescription (ITT]) and not (1701) . The 
series (1751) should be understood as an expansion at large z, possibly asymptotic, of the 
true solution, whose small z expansion is given by ( 1731) . The compatibility of ( 1751) and 
( 1731) is studied more easily for the inverse Laplace transforms. This will be done in the 
next section where we will see that demanding that the two series are compatible fixes 
uniquely the coefficients of both of them. 

4.1.3 Expansion at u = oo and universal scaling function 

By construction, the solution (1661) expands at infinity as 

x W ± 

G ^) = £^T (76) 

n>0 

Comparing the series with the large u asymptotics (JT7J of the physical resolvent 



R phys = -lie + [(1 - iD)G + + (1 + iD)G-) , (77) 

we fix the first two coefficients \(W+ + Wq) = lie and \(Wf + W{) = f{g)/16g 2 . This 
yields a constraint for the expansion coefficients, 

1 = 5>(e) = 5>W a „( e ) , (78) 
and the expression of the universal scaling function f(g) in terms of c n : 

= - + 4j]eH- 1 na n (e). (79) 
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4.1.4 The leading order in the PW/GM limit 

It follows from the scaling (j74p that the solution at the leading order is given by the 
n = term of the series 

G ± (s)=2te (^{) ±4 • (80) 

The constraint (1781) gives Co(0) = 1 and the universal scaling is given by the n = term 
in £f9]): 

f(g) = - e =4g. (81) 

Written for the resolvent and in terms of the variable x(u), the leading order solution 
is 

fl«=o(«) = -2ie | 1 - -j=L= + i-^t=^ ] . (82) 
The density a(u), related to the resolvent by (fl5l) . agrees with the AABEK solution 



4.2 Inverse Laplace transform of the solution 



The relation (1601) involves the shift operator and therefore looks simpler for the Fourier 
transformed quantities. However, in order to be able to exploit the analytic properties 
of the general solution we perform instead an inverse Laplace transformation. Since the 
functions g± and G± are analytic for > 0, we can define the Laplace transformation 
and its inverse 

/(*) = ^ d£ e~ zl f(£) f(£) = -L / dz e zi f(z) . (83) 
Jo z ' 711 Jm.+o 

Similarly we can define the inverse Laplace transformation for the variable z having as 
the origin the left branch point. 

Introduce, similarly to (1641) . the linear combinations 



g± = r+ =F ir- . (84) 
Then from ( IBTl) it follows that the functions g±{z) are related to G± to by 

9± = ^l{D-l)G ± , (85) 

where D = e ldz is the shift operator defined in (1221) . For the inverse Laplace images g± 
and G± this relation takes the form 

- / /a V2sin(|) ~ 
sm § ± | 
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Our aim is to use the relation ( 1861) to investigate the compatibility of the general 
solution (|66|) with the expansion (1731) at z = 0, which in the £-space becomes expansion 
at £ — > oo: 

g ± (£) = r^^ + ^r- 1 . (87) 

n>0 n>0 

It follows from the analytic properties of the resolvents in the rapidity space that, in the 
NFS limit, g±(£) are analytic everywhere except for the negative real axis, while G±{£) 
are analytic everywhere on the positive real axis. We sketch the proof in Appendix A. 
The explicit expression for the inverse Laplace transform of (1751) is a series of confluent 
hypergeometric functions of the first kind 

G ± {£) = ±2i e^a n (e) (n + ±) 1 F 1 (l =p \ T n; 2; -l/e) . (88) 

The PW/GM corresponds to keeping ( = £/e finite when e — > while the NFS 
regime is obtained when keeping £ fixed. In the NFS limit we expand in e with £ fixed. 
Therefore, in order to compare with ( 1H7I) . we are going to use the asymptotic expansion 
the limit £/e — > oo, where the solution has an essential singularity: 

tFt (a; 2; -l/e) oc (£/e)- a 2 F (a, a - 1; ; e/£) /T{2 - a) 

+ e- l ' e {-£/e) a ~ 2 2 F (1 - a, 2 - a; ; e/£) /r(a) . (89) 

The asymptotic expansion of the inverse Laplace image of the solution in this limit is 
evaluated using by f[H9"j) . As far as we interested in the perturbative solution, we can 
neglect the second exponentially small term in flHUj) and write 



G ± (£) =2iJ2 e'-Me V 7 ; 2 F (l T \ T n, tJ T n; ; e/*) . (90) 

In the leading order in e 

G ± (£) = 2i(e/£) 1 ^ (jr ^ + O(e) ] . (91) 

We see that even in the leading order the resolvents scale in the NFS regime as fractional 
powers of e and are linked to the whole perturbative series in the PW/GM regime. In 
the leading order the sum on the r.h.s. of (1911) contains only non-negative powers of £, 
but in the next orders in e more and more negative powers of £ will appear. 
Now we represent, as in [37], the ratio of the sine functions in (186j) as 



sin(|) _ S±(£) 



sin 



2 ± f) T ± (£) 



(92) 



where S and T represent ratios of Gamma functions: 
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If we rewrite the equation (1861) as 

G ± {t) 1 ~g ± (£) 



(94) 



T ± {£) ^2 S ± {£) ' 

then the l.h.s. is analytic everywhere except the negative real axis, while the r.h.s. is 
analytic everywhere except the positive real axis. As a consequence, neither of the sides 
has poles and the only singularities can be branch points at i = and £ = oo. This 
means, in particular, that the expansion of the r.h.s. as a power series at £ = oo coincides 
with the expansion of the l.h.s. at £ = 0. 

To evaluate the coefficients of the two power series we need to expand S± at £ = +oo 
and T± at £ = 0, 



S ± (£) = ±(£/2n)^ 1 + , (95) 




T±{t) = rJT+T) ' 1 •' ! ' (9G) 

As it should, the series expansion of the l.h.s. of ( 1941) at i = contains exactly the same 
fractional powers as that for the r.h.s. at £ = oo. To get rid of these fractional powers, 
we multiply both sides of (194"]) by (£/e) 1=Fl / 4 and write 

S±^W T| =J2 C±{e)£~\ (97) 

where the coefficients C^(e) should be understood as formal series in e, 

C±(e) = Y,C± p e p . (98) 

p=0 

From ( 1H71) . ( 195]) and the relation ( 1941) we deduce that the coefficients in front of the 
non-negative powers of £ vanish, 

C±( e ) = for n = -1,-2,... . (99) 

Solving these contraints (17H1) and ( 1991) order by order in e one can evaluate recursively the 
Taylor coefficients of the series f }T4"|) . The recurrence procedure is possible because at 
each order in e the sum on the r.h.s. of fl9T|) contains only a finite number of negative 
powers of £. 

We have learned that the general solution of the BES equation in the NFS limit is of 
the form 

G ± {£) = (£/ey 1±l *T ± (£) ^C±(e)r", (100) 

n>0 

g ± (£) = V2 V/ey 1 ^ S ± (£)Y, C±{e) £~\ (101) 

n>0 

with computable coefficient functions given by the formal Taylor series (l98l) . Comparing 
the expansions of ( 1901) and (11001) at each order in e and imposing the condition ( 1781 ) we 
can evaluate both sets of coefficients a„ iP and C^ p . We show below how the procedure 
works for the leading order. 
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4.2.1 The leading order in the NFS limit 

In the leading order in e the series expansion of G± at £ = 0, given by (19~T1) . contains 
only non-negative powers in I. Therefore the sum on the r.h.s. of (l9"7j) contains only the 
term with n = 0, and we have 



From the constraint (1751) . which in the leading order gives ao,o = 1, we evaluate 

r(|) _ r(i) 



For the other coefficients we find 



"±n,o = r(±1 . T n , (104) 



where are the coefficients in the expansion 

7r — 6 log 2 



Tr 



4vr 

7r + 6 log 2 
96fsT - 7tt 2 - 36tt log 2 + 108(log 2) : 



2 " 96tt 2 

96K + 77r 2 - 367rlog2-108(log2) 2 

= 64^ ' (1 ° 5) 

(106) 

4.2.2 The universal scaling function 

There is no difficulty to carry out the procedure for the higher orders in e. The only 
diffference will be that the expansion (J9T|) and therefore the r.h.s. of ( 11021) will contain 
some negative powers of I. We do not go into details because the procedure is technically 
identical as the one formulated in E7J . The lowest orders for a„(e) are: 



a„(e) = 1- |e+..., (107) 



7r-61og2 -96K+(77r-121og2)(7r + 61og2) 

ai(e) = ~^6^ + * e + '--' (1 ° 8) 

7r + 61og(2) -96K - 7vr 2 + 547rlog2 + 216(log2) 2 

a - l(e) = 16vr + 38tf e + '--' (1 ° 9) 

5(96K-77r 2 - 367rlog2 + 108(log2) 2 ) 

Me) = ^ + .-., (HO) 



96ir+77r 2 - 367rlog2-108(log2) 2 , . 

«- 2 (e) = ^ + •••• (HI) 

(112) 
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From here we reproduce the result of [37] for the universal scaling function, 



/(e) = - + 4X> h W(e) 



(113) 



n=l 



1 31og2 K 




(114) 



e 7T 7T" 



5 Conclusion 

We have reformulated the Beisert, Eden and Staudacher equation in terms of a functional 
equation obeyed by the resolvent. A similar approach was attempted, although not fully 
exploited, by Kotikov and Lipatov [27] . As shown recently by Basso, Korchemsky and 
Kotahski [37], in the strong coupling perturbative regime it is possible to find the general 
solution as a linear combination of a set of particular functions. This is possible because, 
in the absence of non-perturbative terms, the resolvent possesses extra symmetries. We 
have shown that the "quantization condition" of [37], which allows to fix the coefficients 
of the linear combination order by order in the inverse coupling constant e can be under- 
stood as an analyticity condition on the resolvent. The condition that the resolvent has 
an integrable singularity of the square root type at the points u — ±1, together with the 
conditions on the behavior at infinity of the resolvent are sufficient to fix the solution 
recursively order by order in e. 

Although we have not explicitly investigated the non-perturbative correction, their 
source is clearly identified at the level of the resolvent. This object possesses a sequence 
of self-repeating cuts situated at a distance 2e of one another. When e — > 0, the cuts 
condense and we are left with a single cut plus a non-perturbative term. We leave the 
investigation of the functional BES equation for a future work. 

One of the points of technical importance in the work of [37] and in our work was 
to transform the BES equation into a set of two equations with the so-called undressed 
kernel appearing linearly. In order to perform this transformation we are led to introduce 
an auxiliary density. It would be interesting to know whether such a linearization is 
possible for the general Bethe ansatz equations for M = 4 SYM, and if the auxiliary 
density can be given a physical meaning. Suggestions about the possibility that the 
dressed kernel originates from the elimination of an auxiliary set of Bethe roots have 
been made in [40j |4T] . 

It would be interesting to check if the same method can be applied for the integral 
equations corresponding to other sectors of the M = 4 gauge theory, as well as for the 
other limits in the sl(2) sector. In particular, it would be interesting to try to reproduce 
the new universal scaling function predicted in [16] and [12] and computed by Roiban 
and Tseytlin [43] in the so-called slow long string limit. 
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A Analytical structure of g± and G± 

The inverse Laplace transform f(£), originally defined by (1831) for £ > 0, can be ana- 
lytically continued for complex values of i by rotating the integration contour so that 
asymtotically %t(zi) = at large z. 



X X 

X X 

X X 

X X 

X X 

X X 

-oo ■ -oo • 



X 
X 
X 
X 
X 



Fig. 4: Left: Riemann surface for g±(z). Right: Riemann surface for G±{z). Dots denote 
branch points on the physical sheet, crosses denote the positions of the branch points on the 

lower sheets. 



-oo 




Fig. 5: Deformation of the integration contour for the inverse Laplace transform of G±. 



The functions g±(z) have one second-order branch point at z = on the physical 
sheet and an infinite sequence of equidistant branch points z G iL on the second sheet 
(Fig. 4, left). We assume that we are in the NFS limit in which the left endpoints 
of the branch cuts are sent to — oo. We analytically continue g±{£) beginning from £ 
real positive by changing the phase of £ clockwise. The contour of integration in the 
definition of the inverse Laplace transform will correspondingly rotate counterclockwise. 
For £ G — zR + , the integration contour will lie along the real axis and above the cut of 
g±(z). We can further decrease the phase of £ by rotationg the contour so that half of 
it passes in the lower sheets of the Riemann surface. The procedure can be continued 
without encountering any singularity until the phase of £ is rotated by n 

£ _> e in £ = -£ . (115) 
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At this point, the integration contour goes again on the imaginary axis, with the lower 
half now approaching a sequence of branch points on the lower sheets. Since the contour 
cannot be moved further, we deduce that g±(£) have singularities on the negative real 
axis. 

The functions G±(z) have a sequence of branch points on the negative imaginary 
axis on the physical sheet (Fig. 4, right). Since G± do not decrease sufficiently fast at 
infinity, the inverse Laplace transform does not exist for I > 0. However if we rotate 
slightly the contour counterclockwise, as is shown in Fig. 5, the integral (|83p starts to 
converge. In particular, it is well defined on the negative real axis, when I — > e l7T £. After 
rotating the contour by angle ir, half of it passes on the second sheet, where there are 
no branch points below the real axis. Therefore we can continue rotating the integration 
contour until i — > e 2l7T £, when we encounter the branch points on the second sheet, which 
are on the positive imaginary axis. The inverse Laplace transform G±(£) is therefore well 
defined for £ < and the singularities only occur for I on the positive real axis. 

B Relation with the BKK conventions 

The functions used in the present paper and those used in the paper by Basso, Korchem- 
sky and Kotahski [37] are related as follows: 

oo 

X hcic (u) =-iJ dte Uu Y BKK (t) , 
o 

with 



-^here(w) 




r+(u) 
r_(u) 


7-(0 
7+(0 


r+(u) - 2ie 

r_(«) 


£r-(t) 
-|r + (t) 
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